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Abstract. The present article is the first from a series of three papers which will be 
published in this journal, and which provides a brief presentation of the book, 
“Microwave - numerical solutions”, written by the same author. The main objective of 
the present article is to describe the development of rigorous mathematical models for the 
optimal design of microwave circuits which use microstrip lines. The models will use real 
input data regarding the propagation characteristics of the electromagnetic field and the 
configuration of all existing wave modes in the line and will provide useful results for 
customer specific microwave applications. 
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I. Introduction 


The book “Microwave-Numerical solutions" written by the same author aims to 
make an essential contribution to the elucidation of many of the challenges of the 
electromagnetic field in the microstrip line. In addition, the book highlights a suite 
of interactive program packages designed to study the dynamic behaviours of the 
electromagnetic field using S-parameters in shielded microstrip line and related 
microwave circuits. 

The present article, the first in a series of three papers which will be published in 
the present and the next volumes of the “Annals of Academy of Romanian 
Scientists” is briefly aiming the mathematical description of the electrodynamic 
method applied to the shielded microstrip line. 

In a series of specialized works, [1]+[7], the analysis and calculation of the 
parameters of the microstrip lines is carried out under the assumption of the quasi- 
static approximation, which assumes that the fundamental mode of the 
propagation wave can be approximated with the transverse electromagnetic mode 
(TEM ). Such approach allows to obtain satisfactory results only for the values of 
the longest wavelengths in the microwave frequency range, when the wavelength 
considerably exceeds the transverse dimensions of the line. Recent achievements 
in the field of microwaves require, however, the operation of microstrip lines at 
much higher frequencies, reaching hundreds of GHz, as well as the use of 
substrates with high relative permittivity. As the working frequency increases, the 
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quasi-static analysis of the microstrip line produces increasing errors. This 
phenomenon is the consequence of the dispersive character of the microstrip line 
(parameters vary according to frequency) and the existence of higher order wave 
modes in the line. Since the microstrip line is a non-homogeneous structure, 
containing air and dielectric domains with different properties, the propagation 
mode is a hybrid one and cannot be associated with the TEM propagation mode. 
The study of the behaviours of the electromagnetic field in the shielded microstrip 
line assumes the satisfaction of the following objectives, presented in detail in the 
content of the book: 


1) to deal with the real nature of the hybrid propagation modes and determine the 
higher order hybrid propagation modes and allow to obtain information about the 
dispersion characteristics of the line parameters; 


2) to consider the microstrip line placed inside a metal box, and the conditions 
generated by electrical shielding effects; 


3) to consider, for practical reasons, the fact that the dimensions of the shielding 
box must be much larger, compared to the thickness of the dielectric domain in 
the substrate and the width of the metal strip placed between the air and dielectric 
domains; 


4) to use a sufficiently general method to allow obtaining some general solutions, 
which can be extended to microstrip structures with physical-geometric non- 
homogeneities of the conductors of more complex lines; 


5) to use the right approximations, so that the accuracy of the calculations is 
limited only by the computing power and related software. 


The titles of the next two papers to be published in the coming issues of the 
Annals, and which will briefly cover the major contributions described in the 
book are as follows: 


- The study of the electromagnetic field from the shielded microstrip line using 
the finite difference method; 


- Matlab software package for calculating parameters of electromagnetic field 
and microwave circuits. 


II. Mathematical background of the model 


The analysed method allows verifying the equations of electromagnetism in the 
entire analysed domain and _ satisfying the conditions imposed on_ the 
electromagnetic field at the separation surface between the air and dielectric 
domains in the composition of the screened microstrip line and in the immediate 
vicinity of the edge of the conductor located between them. 
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The Helmholtz equations 


The formulation of the problem aims to satisfy all the objectives presented in the 
introductory chapter and the mathematical formalization of the adopted physical 
model. 


The starting point is represented by the scalar equations for the axial components 
of the electric and magnetic field known as the Helmholtz equations (figure 1): 


A, E,+k3E,=0, (1a) 
M,H,+k}H,=0, (1b) 


and the solutions are determined in the form of series formed by eigen functions, 
which satisfy, by members, equations (la) and (1b), respectively: 


Es (xy) _ Dos Agm X€m (x) Yesm (y), (2a) 
Hz5 (xy) = Ps Bsm Xlm, (x) Yhgm (y), (2b) 
where: 


- Asm and Bs,, are unknown coefficients; 
- Xe,(x) =cosky,_x and Xh,(x)=sink,.x form a system of eigen 


(orthogonal) functions on the interval | 0, s|: 


i Yesm(y) = sin| kysm(y =< bs)| and Yism(y) = cos[ kysm(y ~— bs)| also 


form a system of eigen (orthogonal) functions on the interval [o, = 


- Kym = — yon = kz —k2,,; the index 6 differentiates the two domains in 


figure 1. 
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Figure 1. The elementary cell of shielded microstrip line 


The transverse components of the electric field and the magnetic field are 
determined with the relations: 


Exs=-g (BGS + wttoltrs “5-*) (3a) 
Bys=-73 (B oes — wHloltys a2") (3b) 
Hys5=- = (B we — WEES =) (3c) 
Hys— = (8 ones + weot,5) (3d) 


where &€9&,g and [oft,g are the dielectric domain’s permittivity, respectively 
magnetic permeability (vacuum and relative ones); in propagation conditions, in 
the microstrip line without losses, the propagation constant is purely imaginary, 
respectively: y= i £. 


The method of partial domains 


The method of partial domains has gained wider use in solving the most different 
problems of electrodynamics. 


The practical value of the calculation algorithm is determined by its 
characteristics: the speed of convergence of the chosen solution, the achieved 
precision, and its stability, which depend on the type of eigen functions used to 
approximate the solutions of the electromagnetism equations. 


To elaborate the selection criteria of orthogonal function systems, a specific 
situation is analysed. 
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Figure 2. Model adapted to the possibility of using orthogonal function systems at the separation 
surface between domains I and 2. 


The separation boundary between domains | and 2 (figure 2) coincides with the x 
axis, so that, near the edge, where the electromagnetic field has a particular 
behaviour, the origin of the axes (x=0, y=0) was placed, and in the x plane =1 
there is an electric wall. 


The tangential components of the electric and magnetic fields at the separation 
boundary between domains are used, respectively Fz9, Hz, FExq and Hy», 
expressed in the form of rapidly converging series using orthogonal systems of 


polynomials and Chebyshev functions, which were defined by means of the 
relations: 


0 for O0<x<-> 
E, EZ = 4 
= = Ej for ~<2<* 
2 2 
HH. = Hzz = Hz for -~s2<- (5) 
0 for O0<x<-> 
E Es: = 6 
- - E.9 for eyes (6) 
2 2 
Hey = Hyg = Hy for Sars. (7) 


The tangential components of the electric and magnetic fields at the separation 
boundary between the domains are developed in rapidly converging Fourier 
series, in relation to the system of orthogonal functions given by polynomials and 
Chebyshev functions: 


E,.9(x) — peal CrPen (x), (8) 


Ez9(x) = } DirWen(X), (9) 
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Hyo(x) = Yio Fe Pre), (10) 
Hzo(x) = Yeo GeWne (x), (11) 


where C,,, D,,, Fj, and G, were determined in the contents of the book following 
an extensive calculation process. 


Electrodynamic expressions of the components of the electromagnetic field in 
the shielded microstrip line 


The electrodynamic expressions of the longitudinal and transverse components of 
the electric and magnetic fields in the shielded microstrip line, determined with 
the help of relations (2a), (2b), (3a)+(3d) and (8)+(11) become: 


¥ 71 
Ez5(&Y)=A Dim X€m(x) Soe Yn Gn dn; (12) 
Yesm (v1) 
Ak? Yhgm 
Hz5 (x.y) = Ko a an Xh,, (x ) “i. ; Smal ’ (13) 


where: 
c =p" é _ Bem rN d 
mn O —4in=1 mn€n kj n=1 Onn Tho 


Ac, = Cy, Ady = Dy 


a a = 1207 2. 
€o 


and represents the wave impedance, and the relations (3a)+(3d) become: 


P N 
B Y5m (y) 


eee a.d.— 
kz Yes (1) a eal 


Eys(xsy) =—t AY Xen (2) 


Yhgm 
= Am (x) ae Snns| (14) 


=» # Bem Yesm@) 
Ex5(% y) = —iA Yim Xlm() a ca Ye=1 nn + 


Yhgm(y) Cnnd | (15) 


= Yh5m(¥1) 
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Hyp.) = — ie (2) [Bm 9 ns — 
js : koPobrs ici Yhism (1) m 
= Kg Yesm(9) 
Hrsé E+ ( fe)" Ye’, (a mua 1 %mdn | (16) 
5m () 
H x, =: Z Xhy, les = - WI + 

ao ai 2 aad eI a 

+ (2) em(eerattrs ZEN, and (17) 


The constant A can be determined from the norming condition of the eigen 
functions. In this sense, the normalization condition of the average power 
transmitted through the transverse surface of the shielded microstrip line will be 


used: 
aK xj Hy; — Ey; H;.; )dxdy = 


The components of the tangential electric and magnetic fields at the separation 
boundary between domains, respectively Fz9, Hzg, Exo and H,9, are determined 
with the help of polynomial systems and Chebyshev functions constituted in 
rapidly converging series and defined by means of relations (8)+(11), in which 
replaced the infinite Fourier series with finite partial sums, respectively: 


E,9(x) =A pe Cn Pen (x), 
Ezo(x) =A 2s dn Wen (x), 
Hyo(x) = v1 Fy Qix (x), 


Hzo(x) = Yeas GW), 
The tangential components of the magnetic field, Hz, at the boundary of 
separation between domains, defined by the Cartesian product {y=y7/x FE. | 


expressed in the form of rapidly converging series, were determined in the 
contents of the book following an extensive calculation process. 
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III. Results 


It is used a line configuration, frequently used in complex microwave circuits: 
w=1 mm, a=3,5 mm, yi=0,5 mm, y2=2 mm and €,2= 9 (figure 3). The minimum 
frequency for this line configuration, at which the propagation constant is purely 
imaginary, is 38.188 GHz; it means the lossless propagation in the line occurs for 
working frequencies greater than 38.188 GHz. 

The configuration of the transverse components of the electromagnetic field, at 
the working frequency equal to 38.188 GHz is presented in figures 3 and 4. 

The book also presents the configuration of the longitudinal components of the 
electromagnetic field and other frequencies for which lossless propagation occurs. 


15 imag E,[V/m] 


Figure 3. The variation of the longitudinal component of the electric field in the nodes of the 
network 


4 


x10 


45 real H, [A/m] 205 imag H,[A/m] 


Figure 4. The variation of the longitudinal component of the magnetic field in the nodes of the 
network 


IV. Concluding remarks 


The electrodynamic analysis of the electromagnetic field allows the most accurate 
deciphering of the phenomena existing in the shielded microstrip line under the 
conditions of satisfying all the objectives presented in the introductory chapter. 
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Consequently, with the help of the presented method, the configuration of the 
electromagnetic field, their parameters and dispersion characteristics are 
determined, under the conditions of checking the Helmholtz equations in the 
entirety of the analysed domain, but also under the conditions imposed on the 
electromagnetic field at the separation surface between the two domains and in the 
vicinity of the edge the conductor between them. 

The representations of the electromagnetic field in the shielded microstrip line 
confirm the expectations regarding the fact that its major variations are in the 
domain located in the substrate and in the separation zone between the air and 
dielectric domains. The components H;, Ey and E, are symmetric, and the 
components EF, Hy and H; are antisymmetric in relation to the axis Oy in figure 2, 
in which the elementary cell of the screened microstrip line was represented. 
Another conclusion drawn from the electrodynamic analysis is the fact that the 
propagation of electromagnetic energy in the microstrip line does not occur 
without loss at any operating frequency in the entire microwave range. 
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